Low power sessile droplets actuation via modulated surface acoustic waves 
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Low power actuation of sessile droplets is of primary interest for portable or hybrid lab-on-a-chip 
and harmless manipulation of biofluids. In this paper, we show that the acoustic power required 
to move or deform droplets via surface acoustic waves can be substantially reduced through the 
forcing of the drops inertio-capillary modes of vibrations. Indeed, harmonic, superharmonic and 
subharmonic (parametric) excitation of these modes are observed when the high frequency acoustic 
signal (19.5 MHz) is modulated around Rayleigh-Lamb inertio-capillary frequencies. This resonant 
behavior results in larger oscillations and quicker motion of the drops than in the non-modulated 
case. 

PACS numbers: 47.55. D-, 43.25.Nm, 43.25.Qp, 68.35.Ja 
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One of the challenges in droplet microfluidics is to over- 
come surface capillary forces and contact line retention 
forces, which prevent the motion and deformation of the 
drop [1|. Different techniques such as electrowetting [2|, 
optical toolbox [Sj based on thermocapillary forces in- 
duced by a focused laser, or ultrasonic surface acoustic 
waves (SAW) have been developed to perform simple op- 
erations on droplets. In particular, SAW are efficient 
to achieve actuation, atomization, jetting, oscillations, 
or mixing of small quantities of liquid, either lying on 
a solid substrate or entrapped in confined geometries 
y. However, due to nonlinear coupling between ther- 
mal and acoustical mode, SAW can induce quick tem- 
perature increase in fluid samples [5|, l6[. This can be 
detrimental for the manipulation of biofluids (albumin 
coagulates in a few seconds when excited by SAW), or 
for hybrid SPR (Surface Plasmon Resonance) /SAW lab- 
on-a chip |7| (substrate heating causes a shift in SPR 
reflectivity). To enhance the droplet response, particular 
attention has been devoted to the design of the actuators 
[Si] and to the chemical treatment of the surface [9] to re- 
duce hysteresis and modify the contact angle. However 
less effort has been dedicated to the optimization of the 
acoustic signal in relation to the natural frequencies of 
the drop. Rayleigh in 1879 [io| and Lamb in 1932 [il| 
have identified oscillation modes resulting from a compe- 
tition between inertia and surface tension. While these 
have been first described for levitating drops, they can 
also be adapted to sessile drops. In this case, the droplet 
vibration is affected by the wettability of the surface |12 | 
and the pinning of the contact line [13[ . 

These low- frequency oscillations (typically from 10 to 
200 Hz for millimeter-sized drops) are observed when 
a drop lying on a solid substrate is subjected to sinu- 
soidal SAW of much higher frequency (about 20 MHz) 
\lA ll5| . The free surface deformation results from non- 



linear acoustic forces (acoustic radiation pressure and 
acoustic streaming bulk force), while the detailed mecha- 
nism of excitation of these modes has not been elucidated 
yet. The acoustic radiation pressure induces a stress at 
the surface of the drop, while acoustic streaming induces 
internal flow, which can also contribute to the droplet de- 
formation. The respective magnitudes of these two force 
fields depend on the properties of the acoustic field inside 
the drop and thus on the frequency of excitation, the size 
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of the drop, and the acoustic attenuation length 

In this paper, we show that modulations of the acous- 
tic signal around Rayleigh-Lamb characteristic frequency 
and twice this frequency, result respectively in harmonic 
and parametric response. This latter was predicted the- 
oretically by Papoular and Parayre for levitating drops 
[17| but not observed experimentally. This resonant be- 
havior leads to higher amplitude drop oscillations com- 
pared to the non-modulated case at same input acous- 
tic power. We also show that these oscillations promote 
droplets mobility. 
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FIG. 1: (color online) (a) Sketch of the experimental setup, 
(b) Drop undergoing large deformation along the acoustic 
wave refraction angle (j)r. 

SAW are generated at the surface of a 1.05 mm thick 
piezoelectric substrate (X-cut Niobate Lithium LiNbOs) 
by a transducer consisting of interdigitated fingers. Fig. 
dJ These fingers are designed with the following process: 



1) A titanium (Ti) layer of 20 nm and a gold (Au) layer of 
200 nm are successively sputtered on a LiNbOa substrate 

2) The substrate is coated with AZnlof2020 resist, which 
is patterned by conventional photolithography technique 

3) The Au/Ti layers are successively wet-etched by potas- 
sium iodide (KI) and hydrofluoric acid (HF 50%), and 4) 
AZnlof2020 is removed by acetone. The width of the 
fingers and their distance are both equal to 43.75 /im, 
leading to a characteristic frequency of 19.5 MHz, which 
is used as the carrier frequency /c. A periodic sinusoidal 
voltage is applied at this frequency with a high frequency 
generator (IFR 2023A) and amplified with a home made 
amplifier. This carrier signal is modulated by a square 
wave switching between 1 and 0, at frequency fm ^ fc 
The amplitude d of the SAW is measured with a Mach- 
Zender laser interferometer (BMI-SH130). The surface of 
the substrate is treated with hydrophobic coating (mono- 
layer of OTS) leading to advancing and receding contact 
angles of Oa = 108^ and Or = 99^ respectively, measured 
with a Kruss DSAIOO goniometer. A 7.5 fil droplet of 
water is then placed on the substrate (a sensitivity anal- 
ysis of the droplet response according to its volume has 
been conducted in Ref. Il5l ). The droplet dynamics is 
observed via a high speed camera (Photron SA3) and 
recorded at 2000 frames per second. To avoid pollution 
of the surface by impurities and to obtain reproducible 
results, all the experiments have been carried out in our 
laboratory class 1000 clean room. At rest, the drop re- 
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FIG. 2: (color online) Dipolar (degree 1=2) zonal oscillations 
of a drop of 7.5 /j1 subjected to a SAW of carrier frequency 
fc = 19.5 MHz, modulation frequency fm = 52.5 Hz, and 
amplitude d = 1.38 nm. The time elapsed between two suc- 
cessive snapshots is 2 ms. Ah and Aw are respectively the 
longitudinal and lateral amplitude of oscillation. 

mains essentially hemispherical since its radius {R ^ 1.5 
mm) does not exceed the capillary length (/c ~ 2.5mm 
for water at room temperature) and thus surface tension 
overcomes gravity. The acoustic energy transmitted to 
the fiuid induces zonal drop oscillations, of degree 1=2 
(in spherical-harmonics basis), and a tilt of the drop to 
the right, as shown in Fig. [21 This left-right asymme- 
try is due to the asymmetry of the acoustic field, which 
is radiated in the drop according to a refraction angle 
0^ ^ 25^ given by Snell-Descartes law, see Fig. [TJ When 
the drop is sufficiently tilted for the rear and front contact 
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FIG. 3: (color online) Amplitude of vertical oscillations of 
the drop Ah divided by the initial height of the droplet ho as 
a function of the modulation frequency fm for different am- 
plitudes d of the surface acoustic wave. In the green, blue 
and red region, the droplet response is respectively super- 
harmonic, harmonic and subharmonic (compared to the fre- 
quency of modulation) . 



angles to exceed their hysteretic value, one additionally 
observes drop motion. The vertical amplitude Ah and 
frequency fr of oscillation depend on both fm and d^ see 
Fig.[3l According to these parameters, either a harmonic 
response at frequency fr = fm^ ^ subharmonic response 
at fr = \l2fm or a super harmonic response at fr = 2/^ 
is observed (see also accompanying movies in the supple- 
mental material). 

Harmonic region: shift of resonance. In the blue 
area of Fig. [3l a peak of response is obtained when fm 
reaches the inertio-capfflary characteristic frequency /o, 
which can be estimated from Rayleigh formula: fo = 
(87/37rpV)^/^ ^ 89 Hz for the dipolar (/ = 2) oscilla- 
tion of a 7.5/i/ droplet, with 7 the surface tension and V 
the droplet volume. At intermediate power ((i=1.10 nm, 
(i=1.38 nm), the peak is asymmetric, with a skewness di- 
rected to low frequencies and the resonance frequency de- 
creases with the amplitude of oscillation. This response is 
typical of an anharmonic oscillator with softening spring 
W < 0): 



• 2Xx + uj^x + ax'^ 



px^ = Fcos{nt) 



(1) 



where x is the dynamic variable (here the deformation 
of the drop), t the time, A the damping coefficient, 
ujo = 27r/o the angular eigen frequency, a and P two non- 
linearity coefficients, F the amplitude of excitation and 
Q the excitation frequency. Such nonlinear behavior has 
already been reported by Perez et al. [l8[ for levitating 
drops larger than the capillary length and more recently 
by Miyamoto et al. [l9[ for sessile droplets smaller than 



the capillary length, with pinned contact line. These au- 
thors determine the coefficient appearing in Eq. ([T]) from 
experiments and compare the frequency response of the 
drop to theoretical predictions. In our system the oscilla- 
tion damping is due to dissipation in the viscous bound- 
ary layer and in the neighborhood of the contact line [20| . 
The nonlinear response of droplets appears when they 
undergo finite- amplitude deformations [2l|, |22|, leading 
to a shift of the resonance frequency to lower frequen- 
cies as the amplitude of oscillation increases. For sessile 
drop, additional nonlinearity results from the up-down 
asymmetry of boundaries and the presence of a contact 
line. 

Superharmonic region: combination of modes. In 
the green region of Fig. [H low fm and large A/i, the 
drop response is a combination of harmonic and super- 
harmonic modes. Indeed, the drop is pushed by the 
acoustic wave when the signal is on. Then, the drop 
keeps bouncing in the period with no forcing (signal is 
off) and oscillates a whole cycle before the next push. Af- 
ter a transient phase, this synchronization of forced and 
natural bouncing results in large drop oscillations. 
Subharmonic region: parametric resonance. In 
the red region of Fig. [3l the droplet responds at /m/2. 
This subharmonic response appears only above a thresh- 
old: d > 1.42 nm, with a frequency window broaden- 
ing progressively with the amplitude of oscillation Ah. 
This is typical of a so-called Arnold Tongue. Further- 
more, the amplitude of the subharmonic response at fixed 
fm is independent of the amplitude of excitation d. All 
of these properties are characteristic of parametric reso- 
nance. Parametric instability of an oscillator is enabled 
when its characteristic frequency uJq is modulated in time 
near 2uJo [23']. It can be modeled with a Mathieu equa- 
tion: 



X + 2Xx 



■^lii 



Acos{2ujo + e)t]x 



(2) 



with 6 <C cJq. Such Mathieu equation can be obtained 
when an anharmonic oscillator described by Eq. ([T]) is 
excited near 2(jOo: ^ = 2ujo + e [2J]. Indeed, from the 
asymptotic expansion of the variable x in Eq. ^: x = 
xi + ^2, with xi the solution of the harmonic oscillator 
and X2 <C xi, we obtain at second order: 



X2 
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3^ 



cos(2cc;o + e)t - 



(3) 



The parametric excitation of zonal oscillation modes was 
predicted theoretically by Papoular and Parayre [17| for 
levitating drops but not observed experimentally. In- 
deed, parametric resonance appears for anharmonic os- 
cillators above an amplitude threshold Ft = 6Acj^/|a| 
which decreases with the nonlinearity coefficient a. This 
latter is related to the asymmetry of the oscillator stiff- 
ness for prolate and oblate deformation (x < or > 0) 
which is increased by the presence of the substrate. Thus, 



the observation of the parametric response in the present 
experiments could be explained by a larger nonlinearity 
of the system and larger excitations than in the system 
of Papoular and Parayre. 
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FIG. 4: (color online) Evolution of a 7.5 /A droplet excited by 
an acoustic wave of amplitude d = 1.55 nm and modulation 
frequency fm = 100 Hz. The top picture shows the initial 
shape of the drop. The spatio-temporal diagram below is 
obtained by taking the base line of the drop (dashed blue line 
on top picture) and showing its evolution as a function of 
time. The average slope of the front or rear curve delimiting 
the drop gives the velocity of the drop. 

The parametric mode appears after a transient phase, 
as shown in the spatio-temporal diagram on Fig. HJ At 
first, the droplet oscillates at the same frequency as the 
excitation fm = 100 Hz between t = and 0.15 s. Be- 
tween t = 0.15 and t = 0.32 s, the parametric insta- 
bility grows. Finally, the droplet reaches a stationary 
regime of oscillation at half the frequency of modulation 
fr = l/2fm- Between initial harmonic and later para- 
metric response, the amplitude of droplet lateral oscilla- 
tions Aw is increased by a factor of 4.8. This large in- 
crease of Aw comes along with an increase of the droplet 
velocity by a factor of 4, see the rupture of slope at t = 
0.35 s in Fig. H 

In all the previously described experiments, there is 
indeed a correlation between the droplets velocity V and 
their amplitude of oscillation Ah/ho^ see Fig. [5l Indeed, 
the date of the drop velocity divided by the applied acous- 
tic power collapse into a single curve, which increases 
with the amplitude of oscillation up to Ah/ ho = 1, and 
then reaches a plateau. From this figure, we can con- 
clude that (i) the droplet velocity is basically propor- 
tional to the surface acoustic wave power and (ii) that 
droplets oscillations promote their mobility. Indeed, the 
velocity of a drop is determined by the equilibrium be- 
tween the driving forces and the retention forces. The 
motion is due to the contact angle difference between the 
front and rear interfaces and thus, to the asymmetry of 
the drop, see e.g. Ref. |20|. An increase of the acoustic 
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FIG. 5: (color online) Velocity of the drop divided by the 
square of the surface acoustic wave amplitude d^ (correspond- 
ing to the acoustic power radiated into the drop up to a prefac- 
tor) as a function of the longitudinal amplitude of oscillation 
Ah/ ho. Each marker corresponds to a specific amplitude d. 



FIG. 6: (color online) Velocity of the drop V measured with a 
modulated signal divided by a reference velocity Vref obtained 
at same acoustic power when the modulation is turned off, as 
a function of the modulation frequency fm ■ 



wave amplitude d naturally leads to larger droplet defor- 
mations and therefore larger asymmetry and velocities 
of the drop. The signal modulation does not affect the 
acoustic power. Nevertheless, it can affect both the driv- 
ing and retention forces into different ways. First, it is 
important to mention that while some authors, see e.g. 
Ref. |25|, have observed a decrease of the hysteresis of the 
contact line due to its continuous depinning, such varia- 
tions are not observed in the present study and cannot 
explain the increase of drop velocity. Second, as seen pre- 
viously, the droplet static shape is not much affected by 
gravity since the droplet radius is smaller than the capil- 
lary length. However, when the drop is highly stretched 
with a left/right asymmetry component (see e.g. Fig. 
lb), gravity plays a significant role during the retraction 
phase. The drop is not simply pulled back by capillary 
forces, but also the upper part of the liquid drop falls 
vertically by gravity (see accompanying movie "super- 
harmonic" in supplemental material): this strengthens 
even further the left/right asymmetry and should con- 
tribute to the drop motion. Finally, the acoustic energy 
transferred to the drop and the retention force depend re- 
spectively on the contact surface with the substrate and 
the perimeter of the contact line. Thus a stretching of the 
drop contributes to a reduction of the retention force but 
also a reduction of the amount of energy transferred to 
the drop. This might explain the existence of the plateau 
region in Fig. [5l To conclude, we can compare the veloci- 
ties obtained with and without signal modulation at same 
acoustic power (see Fig. [6] and accompanying movies in 
the supplemental material). The drop velocity can be 
increased by a factor 100 thanks to the modulation. 

In this paper, we have shown that inertio-capillary 
modes of oscillations can be excited either directly 



or parametrically by modulating the acoustic signal 
around once or twice the Rayleigh-Lamb characteristic 
frequency. This modulation introduces an original way 
to decrease the acoustic power required to stir inside, to 
stretch or to move a sessile droplet with SAWs. Com- 
pared to the non-modulated case, the minimum acoustic 
power required to move a droplet at a non-zero speed, 
and a speed of 5 mm s~^ are reduced by a factor of 2 
and 3 respectively, while the minimum power required 
to stretch the droplet vertically by a factor Ah/ ho of 0.2 
and 1 are reduced by a factor 5 and 3 respectively. This 
is especially important for the harmless manipulations 
of biofluids which could be damaged by the temperature 
increase due to the dissipation of acoustic energy. It is 
also of primary interest for portable lab-on-a-chip which 
requires low power consumption. 
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